o 
O 

m 



X 



On the possible observational manifestation of 
supernova shock impact on the neutron star 

magnetosphere 



oo 

O ; A.E. Egorov, K.A.Postnov * 

Sternberg Astronomical Institute, Universitetski pr. 13, Moscow 119992, Russia 



Abstract 



Impact of supernova explosion on the neutron star magnetosphere in a massive binary 
system is considered. The supernova shock striking the NS magnetosphere filled with 
plasma can lead to the formation of a magnetospheric tail with significant magnetic energy. 
The magnetic field reconnection in the current sheet formed can convert the magnetic 
energy stored in the tail into kinetic energy of accelerated charged particles. Plasma 
Q , instabilities excited by beams of relativistic particles can lead to the formation of a short 

5— i ■ pulse of coherent radio emission with parameters similar to those of the observed bright 

' extragalactic millisecond radio burst (Lorimer et al. 2007). 



Subject headings: stars: neutron — radio: bursts — pulsars: general 



'. 1 Introduction 

| Supernova (SN) explosion in a binary system with neutron star (NS) is not a rare event. As 

{SJ ■ the binary system can be disrupted during the second SN explosion, the lower limit to the rate 

of such events can be estimated, for example, from the galactic binary NS formation rate. This 
rate as inferred from binary pulsar statistics and population synthesis calculations is of order 
OO ! °f 10~ 4 — 10~ 5 per year (Postnov & Yungelson 2006), i.e. a few - ten events per cubic Gpc per 

day. 

The impact of the SN shock on NS magnetosphere can lead to spectacular observational 
manifestations. For example, Istomin & Romberg (2002) argue that it can result in a narrow- 
collimated jet of hard radiation which can be observed as a gamma-ray burst (GRB). Their 
arguments are based on the cumulative effect of the SN shock acting on the NS magnetosphere 
resulting in the formation of a long magnetospheric tail. The magnetic field reconnection in 
the tail then gives rise to particle acceleration. Synchrotron emission of these particles in the 
magnetic field can produce hard X-ray photons observed as a gamma-ray burst. Conversion of 
magnetic energy in the tail into hard X-rays is a complicated physical problem. It is still to be 
seen if this mechanism can explain the gamma-ray burst phenomenon. However the generation 
of relativistic particles during magnetic field reconnection seems to be unavoidable (see recent 
numerical calculations in (Zenitani & Hoshino 2007)), and a large variety of plasma phenomena 
are expected to occur during interaction of these particles with the ambient medium. 

Here we suggest that plasma instabilities in the process of the magnetic field reconnection 
can be responsible for the formation of a powerful short non-thermal radio burst. A very short 
radio burst with peculiar properties was recently discovered by Lorimer et al. (2007). The 
radio flux of the burst at frequency 1.4 GHz was 30 ± 10 Jy, the duration of the event was 
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estimated to be shorter than 5 msec and its dispersion measure was 375 cm _3 pc suggesting 
the possible extragalactic nature of the source at a distance of ~ 1 Gpc. No host galaxy up 
to the 18th i3-magnitude has been found implying a distance limit of >600 Mpc for a Milky 
Way-like host. Such a distance implies a radio energy release in the burst of ~ 10 40 erg with a 
brightness temperature of ~ 10 34 K for isotropic radiation. The total rate of such bursts was 
estimated roughly to be around 90 d _1 Gpc -3 , which is much lower than the core-collapse SN 
rate ~ 1000 d -1 Gpc~ 3 , but well in excess of the GRB rate ~ 4 d -1 Gpc -3 and the expected 
rate of binary neutron star coalescences ~ 2 d -1 Gpc -3 . This estimate, of course, is very 
uncertain due to a small statistics. 

We should note here that properties of the short radio burst reported in Lorimer et al. 
(2007) are in fact very similar to those of pulsar radio emission: the coherence, an extremely 
high brightness temperature at frequency v rts 1.4 GHz, a steep spectral shape (~ v~ A ). Lorimer 
et al. (2007) discuss several possible known sources of the strong millisecond radio bursts, 
including rotating radio transients (RRATs) and giant pulses from radio pulsars, but none of 
them appear to be energetic enough. Based on statistical arguments only, Popov & Postnov 
(2007) proposed that this burst could be produced by a hyperflare of an extragalactic magnetar. 
Also it was estimated in (Popov & Postnov 2007) that assuming the ratio of radio emission 
power to the magneto-rotational losses as observed in giant radio pulses in Crab with a peak flux 
of 1000 Jy, a simple scaling leads to the conclusion that magnetar magnetic fields (B ~ 10 15 G) 
and millisecond rotation periods are required. However, magnetars are observed to rotate slowly 
P = 1 — 10 s due to rapid spin-down. In the ms magnetar model (Usov 1992, Usov 1994) the 
characteristic spin-down time is several seconds, so if radio emission were produced in a newly 
born magnetar by some mechanism, we would observe a series of ms pulses with decreasing 
period, not a single bright burst. 

In contrast, the magnetic field reconnection during relaxation of NS magnetosphere after 
impact of SN shock naturally leads to the formation of relativistic particle streams. Plasma 
instabilities in such flows stimulate conversion of kinetic energy of particles to the energy of 
transverse electromagnetic waves, which can be observed as the short powerfull burst of non- 
thermal radio emission. Lyubarsky (2008b) has shown that such radio bursts can propagate 
through the stellar wind without being destroyed by the induced Compton and Raman scatter- 
ing provided that they originate at distances > 10 15 cm from the exploding massive star and 
their source is moving relativistically. While in massive binary systems such distances from 
the SN progenitor may be natural, the relativistic bulk motion of the emitting region requires 
justification. In our model emitting region moves relativistically due to acceleration of electrons 
during non-stationary magnetic reconnection. 



2 Shock interaction with NS magnetosphere 

Following Istomin & Romberg (2002), we consider typical SN shock parameters at the free 
expansion stage after the explosion (Imshennik & Nadyozhin 1989): the total kinetic energy 
Ekin ~ 10 47 erg, the velocity u « 4 • 10 9 cm/s, the density p w 10 -8 g/cm 3 , and the shock 
width h w 10 9 cm. 

Assume that the NS was a pulsar just before the explosion, i.e. relativistic plasma was 
generated near its surface. A young NS can shine as radio pulsar for about 10 7 yrs. In a massive 
binary, the probability for a SN explosion during this stage is not small, especially if the original 
binary components had close masses. The characteristic size of the NS magnetosphere (ignoring 
possible asymmetry due to a strong stellar wind of the secondary component) is about the light 
cylinder size Ri ~ c/fl = cP/2n. The NS magnetosphere must be filled with pair plasma with 
the number density determined by the Goldreich- Julian value (Goldreich & Julian 1969): 

\pG.i\k \flB\k 
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where fc = 10 4 — 10 5 is the particle multiplication factor, e is the elementary charge. 

The shock kinetic energy geometrically intercepted by the NS magnetosphere can be esti- 
mated as 

nR 2 , „ M ( A ~ 2 / a \ -2 

lOrad/s, 

where a is the distance between the stars. According to Istomin & Romberg (2002), such 
shock perturbation can lead to the formation of a magnetospheric tail with substantial magnetic 
energy. The tail is unstable: magnetic reconnection occurs and charged particles are accelerated 
to high energies. We suggest that the disruption of the magnetospheric tail can lead to the 
formation of two colliding relativistic electronic and positronic beams. Plasma instabilities 
developed in the collision generate plasma waves which further can be converted into transverse 
electromagnetic waves, giving rise to a narrow-collimated beam of outgoing non-thermal radio 
emission. The situation here is similar to the most plausible mechanism of radio emission 
generation in pulsars (Lyubarsky 1992, 2008a). 

Istomin & Romberg (2002) argue that SN shock passing through the NS magnetosphere 
compresses its head-on part and forms a long narrow magnetospheric tail like in the terrestrial 
magnetosphere blown by solar wind. In such a magnetospheric tail conditions of the magnetic 
field freezing into the shock plasma and the magnetic flux conservation in the tail (see Istomin 
& Romberg (2002) for more detail) imply that the magnetic field strength in the tail and the 
tail's diameter at distance I from the NS are 
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respectively. Here 
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dt = 2r* ( ^ ) , (4) 
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B* = (Anpu 2 ) 1 ' 2 « 10 6 5 - " , G 

V H ' ^lO-Vcm 3 ; V4-10 9 cm/s/ 

is the magnetic field at the Alfvenic radius r* determined as a distance where the magnetic field 
pressure is balanced by the shock ram pressure. For the dipole magnetic field and the assumed 
shock parameters we find 

- " R {i ) " -* (j£0 " 6 (loL) (to^g ) " - < s > 

where B is the field strength at the NS surface, R is the NS radius. So the magnetic energy 
storage in the narrow (dt <C L ~ h, where L is the tail length), magnetized (B* rts 10 6 G) 
magnetospheric tail will be 

« - / = r^'^ - ^ die) 2 (^r (i^) 3/2 - O 

Clearly, such a tail can be unstable with respect to magnetic reconnection in the current sheet. 
The symmetry of the problem suggests the plane geometry of the sheet in the orbital plane. 
During the magnetic reconnection, the electric field component transverse to the tail axis is 
generated in the plane sheet. This field can accelerate plasma electrons with a broad energy 
distribution up to high Lorentz factors (see e.g. Zenitani & Hoshino (2007)). Importantly, 
some amount of electron-positron plasma is initially expected to be present in the tail with 
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the number density below the Goldreich- Julian value ([T]), since during the formation of the 
magnetospheric tail electron-positron plasma will outflow along the magnetic lines that have 
been opened up by the shock. The upper limit on the number density of pair plasma in the tail 
is 

n pairs < noAD « 10 7 (^) (^) (A) 3 cm -a. (7) 

At the same time, some amount of electron-ion plasma can penetrate from the shock into the 
tail during its formation. The number density of electron-ion plasma can significantly exceeds 
that of pairs. 



3 Generation of radio emission 

At the beginning of the magnetic reconnection, the tail is filled with plasma with admixture of 
pairs. During the reconnection, several X-points can be formed in the long tail (tearing insta- 
bility) and a strong electric field directed transversely to the tail axis appears and accelerate 
charges. Electrons and positrons are accelerated in opposite directions so a low-density positron 
beam can appear propagating against relativistic electron plasma. The situation becomes sim- 
ilar to the pulsar wind flowing along the open field lines of the NS magnetosphere (Lyubarsky 
1992), where tge electric field appears due to leak of Goldreich- Julian charge density (quasi- 
stationary violation of condition of full screening of electric field). A fraction of the kinetic 
energy of accelerated particles can be transformed by plasma instabilities into plasma oscilla- 
tions mostly at the comoving plasma frequency. These oscillations can then be non-linearly 
converted into transverse electromagnetic waves observed as pulsar radio emission. 

In our case, the duration of the radio pulse produced during magnetic reconnection (fast 
reconnection with Alfven velocity) can be roughly estimated from the transverse size of the 
emitting tail: t rs r*/v a ~ r*/c ~ 3 ms. This time is about as observed in the reported 
extragalactic radio burst (Lorimer et al. 2007). 

The subsequent picture can be as follows. During the magnetic field topology changing, 
a stream of relativistic electronic plasma moving transversely to the tail axis with some ef- 
fective mean Lorentz-factor T can emerge. The reverse flux of positronic plasma accelerated 
by the same electric field in the reconnecting layer moves with about the same Lorentz-factor 
against the electrons. Due to various instabilities (beam instability etc.) the kinetic energy of 
the positronic beam dissipates in the electronic plasma by exciting plasma oscillations. Such 
energy conversion is effective in the case of small number density of positrons in comparison 
with electrons (n + <^ n~), which is expected in our case. These oscillations are non-linearly 
converted into transversal electromagnetic waves which can be observed as a short radio burst 
for an appropriate orientation of the tail and the observer. Relativistic aberration confines 
radio emission within a narrow cone with opening angle 9 ~ T -1 . The oscillations occur at the 
comoving plasma frequency 

/ 47re 2 n~ 

The conversion rate of plasma oscillations into transversal waves can be expressed as (e.g. 
Lyubarsky (1992)) 

^ = aW t W h (9) 

where 

&=~~ 2 J (10) 

n mc z 

Wt and Wi are the energy densities of the transversal and longitudinal waves, respectively, in 
the comoving plasma frame. The observed frequency is Doppler boosted (emitting electrons 
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move towards the observer with Lorentz-factor T): 

v = 2TD = 2T^- = 2Te\l — , (11) 

Z7T ' 

where rie is the electron number density at the moment of emission. It can be slightly less than 
the initial number density in the tail since during the magnetic field reconnection the volume 
of emitting current sheet can widen due to some instabilities (Zenitani & Hoshino (2007)). We 
shall adopt rie — an~ , where a < 1 is the ratio of electron number densities in the current 
layer and in the tail. The total energy of generated radio emission E can be estimated as 
the kinetic energy of positrons enhanced by factor (2I 1 ) 2 . One power of factor (2r) appears 
due to the Doppler effect, and the second power is due to transformation of the kinetic energy 
of positrons to the frame co-moving with electronic plasma. Below we shall not distinguish 
between Lorentz-factors of accelerated electrons and positrons. Then 

E = (TN + mc 2 (2T) 2 w 4(T 3 n+7:r* 2 hmc 2 , (12) 

where the factor C < 1 takes into account the efficiency of the kinetic energy conversion into 
electromagnetic radiation, the efficiency of particle acceleration etc. The N + is the total number 
of positrons, ft is the length of the tail. The observed energy flux is then 

E _ 16(r 5 n+r* 2 hmc 2 

TA iy7 r(O.50r) 2 ~ rA^ ' (13) 

where = 1/T is the relativistic aberration angle, A„ is the frequency band of the receiver, r 
is the distance to the source. 

Taking the millisecond radio burst parameters i/ w 1.4 GHz, S ~ 30 Jy, r ~ 500 Mpc, 
A„ w 300 MHz (Lorimer et al. 2007) and assuming n+ = (3n~,f3 < l,r ~ r*/c from Eq. 
(ll)-JH we find 
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The obtained value of the Lorentz-factor does not contradict to constraints imposed by the 
induced scattering in the source plasma (Lyubarsky 2008b). The positronic plasma number 
density and the total energy of radio emission accord with Eqs. (7) and (6), respectively. 

The spectrum of the observed ms radio burst is dS/dv <~ z/~ 4 . As we mention above, the 
generation mechanism of the burst is analogous to that of pulsar radio emission. Consequently, 
it is natural that observed spectrum is similar to pulsar radio emission (at the high frequency 
range, particularly the Crab pulsar has very similar spectrum). Many physical factors can affect 
the spectral shape: the wide spectrum of plasma waves (Aoj ~ uj), the wide energy spectrum 
of accelerated particles, plasma inhomogeneities, plasma refraction near the source, extinction 
in the parent galaxy, etc. The observed spectrum can be partially explained by assuming the 
power-law energy distribution of accelerated particles dn~ /dq ~ 7~ p . The contribution from 
electrons with Lorentz-factors [7,7 + dj] into the observed emission in linear approximation 
can be written as 

dS ~ -fdn- ~ 7 1 - p rf 7 . (18) 

The additional Lorentz-factor is due to the Doppler effect. Then taking into account (11) and 
(18) we obtain 

dS dS ,dv , , 

= / ^ ^ ' ^ 1/ * 

dv c?7 c?7 

This spectrum is in agreement with observations provided that p = 4 — 6, which generally does 
not contradict to the energy spectrum of accelerated non-thermal particles found in numerical 
simulations of relativistic magnetic reconnection (Zenitani & Hoshino (2007)). 

There could be a problem with the generated radio burst passing through the shock. How- 
ever, the shock can be inhomogeneous and there can be regions with significantly lower densities 
(holes) transparent for radio emission. The upper limit to the allowed density p' can be esti- 
mated from the condition that the plasma frequency of the shock matter (fully ionized) be less 
than the frequency of radio waves: 

, Trm p m e v 2 14 3 

P ^ m P n e = — 2 ~ 10 g/cm < p, 

where the n' e is the electron number density corresponding to the maximum allowed plasma 
frequency. Indeed, observations of SN1987a suggested strong inhomogeneities in the expanding 
envelope allowing the early appearance of X-ray emission from the SN remnant (Grebenev & 
Syunyaev 1987). The radio emission can also be seen through inhomogeneities, but the observed 
duration and flux should be modified by effects of passing through such an envelope. 



4 Hard radiation generation 

Let us roughly estimate properties of the possible hard synchrotron emission which can ac- 
company the magnetic reconnection in the magnetospheric tail. The characteristic energy of 
synchrophotons is 

V10 6 G / V2000 J 

which falls into the soft gamma-ray energy range. The upper limits to the total energy of this 
radiation E 1 and the flux S 1 can be estimated from Eqs. (3), (14), (16) with account for the 
synchrotron intensity from one particle I — (2/3)(e 4 Bf)/(m 2 c 3 )T 2 : 



hw~, ~ h — -T z w 40 I ) I ) keV , 

mc 



E 7 = IN- = ^J^n-nr* 2 h « 4 • 10 30 ( -A_V ( — 
7 3m 2 c 3 V 106G / V 2000 



3 • 10 4 cm- 3 / V10 8 cm / V lo9cm / V3ms/ 



erg, (19) 
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-^V (20) 
s • cirr 

Apparently, the total energy of hard non-thermal radiation is similar to that of radio emission 
(17), but the flux falls short of the sensitivity of modern detectors. So no accompanying 
"gamma-ray burst" is expected, which agrees with the absence of any gamma-ray event at the 
time of the observed radio burst (Lorimer et al. 2007). Clearly, should this happen in our 
Galaxy or in the Local Group, both radio and gamma-ray bursts can be observed, probably 
associated with an emerging supernova. 



5 Discussion 

We have shown that a bright millisecond burst of non-thermal radio emission with properties 
similar to those observed by Lorimer et al. (2007) can be explained by the supernova shock 
interaction with NS magnetosphere in a massive binary system. The shock passing through 
the magnetosphere leads to the formation of a long magnetospheric tail (Istomin & Romberg 
2002). We conjecture that the magnetic energy stored in the tail can be transformed into the 
kinetic energy of non-thermal particles due to magnetic reconnection. If there is some fraction of 
positrons (which seems reasonable if NS before the SN explosion was at the pulsar stage), plasma 
instabilities can transform the kinetic energy of the less dense positronic beam into the plasma 
oscillations in the electron beam which can then be non-linearly converted into transversal 
electromagnetic waves, like in pulsar magnetospheres (Lyubarsky 1992), thereby producing a 
short pulse of non-thermal radio emission. The magnetic fields of young NSs can vary in a wide 
range from 10 11 to 10 15 G, leading to a broad range of magnetospheric sizes, densities of pair 
plasma filling the magnetosphere, etc. The situation is even more diverse in massive binary 
systems with strong stellar winds from presupernova stars. In addition, the density and velocity 
of SN shocks can also vary significantly. So all basic parameters of our model (r*, n+, t, V and 
n~) can be very different, and so should be the properties (frequency, flux, duration) of possible 
radio transients accompanying the magnetosphere topology restructuring during reconnection. 

Magnetic reconnection can also occur in the region with toroidal magnetic field beyond the 
light cylinder of a pulsar. In this case the current sheet is located in the rotational equatorial 
plane of the pulsar. The magnetic field strength at the light cylinder Bi can be significant: 

^(l)^-^^)^)^. 

Hence there is a substantial magnetic energy storage cb ~ (Bf /8ir)Rf ~ 10 35 erg outside the 
magnetosphere (taking into account that toroidal magnetic field decreases slowly as 1/r). A 
stationary reconnection of the toroidal field was suggested by Lyubarsky (1996) to explain hard 
emission of pulsars. In our scenario, the passing shock compresses the NS magnetosphere and 
can initiate a non-stationary large-scale reconnection of the toroidal magnetic field beyond the 
light cylinder. Estimates similar to those given in Section 3 can explain the generation of short 
radio burst in two-dimensional current sheet approximation. Its properties can be similar to 
those observed by Lorimer et al. (2007) assuming model parameters similar to those given by 
Eqs. (14)-(17) (provided that positrons, which excite plasma waves, are collected from a large 
volume ~ Rf). Such parameters do not contradict all existing constraints, so this approach 
can be considered as an alternative model for the millisecond radio burst. In this case it is not 
needed to explain the transparency of the shock to the emergent radio emission since the radio 
pulse can be produced before the SN shock has passed the magnetosphere and screened the 
source. 
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The search for other short radio transients is under way (Lorimer et al. 2007), and larger 
statistics can be used in future to test the proposed model. 
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